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Abstract 

A (/3-exponential distribution is a generalization of an exponential distribution 
associated to functions ip in an appropriate class, and the space of (/j-exponential 
distributions has a dually flat structure. We study features of the space of ip- 
exponential distributions, such as the convexity in Wasserstein geometry and the 
stability under an evolution equation. From this study, we provide the new charac- 
terizations to the space of Gaussian measures and the space of g-Gaussian measures. 



1 Introduction 

A Gaussian measure is an exponential distribution on with mean and covariance matrix 
parameters. The space of Gaussian measures has good behaviors such as the flatness in a 
dual structure which is a Riemannian metric with a pair of connections being orthogonal 
to each other, the convexity in Wasserstein geometry which is a metric geometry on the 
space of probability measures, and the stability under the linear evolution equation called 
the Fokker-Planck equation. As for the validity of these behaviors, for instance, the dually 
fiat structure on the space of Gaussian measures leads the Cramer-Rao lower bound for 
any unbiased estimator. 

Recently, Naudts [4j introduced the space of probability measures, which is a gener- 
alization of the space of Gaussian measures equipped with a dually fiat structure, and 
generalized the Cramer-Rao lower bound, where each probability measure is determined 
by the inverse function of the function t t— )■ l/ip{s)ds on (0, oo) for an increasing, posi- 
tive, continuous function (p on (0, oo). The probability measure is called a ip- exponential 
distribution and coincides with a usual exponential distribution when (/^(s) = s. The 
(/3-exponential distribution on has a finite second moment if (p has a certain limiting 
behavior at and oo, which corresponds to the instances of functions if given by ^{s) = s'^ 
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with q G (0, + 4)/((i + 2)). When ip{s) = s'^ except for q = I, the y^-exponential dis- 
tribution with mean and covariance matrix parameters is called a q- Gaussian measure, 
which is a power-law distribution. 

The space of g-Gaussian measures behaves well due to the convenience of power- 
law distributions. Indeed, the space of g-Gaussian measures is convex in Wasserstein 
geometry and stable under the nonlinear evolution equation of porous medium type (see O 
Proposition 5] and [7, Theorem A]). We prove that the convexity and the stability are 
features only of the space of Gaussian measures and the space of g-Gaussian measures. 
In other words, the condition ip{s) = s'^ with q G (0, {d + 4:)/{d + 2)) is a necessary and 
sufficient condition for the space of 99-exponential distributions on M'^ with mean and 
covariance matrix parameters to have each of the two features. 

The paper is organized as follows. In the next section, we give definitions and prop- 
erties about Wasserstein geometry and the yj-exponential function. Section 3 concerns a 
condition for (y9-exponential distributions to have a finite second moment ( Prop osition 13 . 2 p . 
The two spaces of (^^-exponential distributions with mean and covariance matrix param- 
eters are introduced in Section 4, one of which has a dually fiat structure and the other 
is convex in Wasserstein geometry (Proposition 14. 3p . We discuss when the two spaces 
coincide with each other (Theorem 14. ip . In Section 5, we consider the stability of the two 
spaces under a certain evolution equation (Theorem 15. II and Corollary 15.21 ). 

Acknowledgements. The author would like to thank Shin-ichi Ohta for many valuable 
advice and discussions. She is also grateful to Sumio Yamada for helpful comments. She 
is partially supported by JSPS-IHES (EPDI) fellowship. 

2 Preliminaries 

2.1 Wasserstein geometry 

Let us briefiy recall properties of (L^-)Wasserstein geometry over Mf^. See [H |9] and 
references therein for further information. 

In this paper, any measure is always supposed to be a Borel measure. A measure n 
on Mf^ is said to be absolutely continuous with respect to the Lebesgue measure on 

if there exists a measurable function p on M.'^ such that /i = pC^. For a measurable 
map F on and a measure p on M.'^, the push-forward measure F^p of p thorough F is 
a measure on M'^ defined by Fij^p[B] := p[F^^{B)] for all Borel sets B C M'^. Given two 
measures p = pC^, v = aC^ and a differentiable map F on M'', z/ = F^p is equivalent to 
p = a{F)det{dF) /i-almost everywhere, where dF is the total differential of F. For any 
probability measures p and u on W^, a coupling it of p and u is a probability measure 
on X R'^ with marginals p and u. In other words, we have 7i[B x M*^] = p[B] and 
B] = ulB] for all Borel sets B C M.'^. We denote by V2 the set of probability 
measures p on M.'^ with finite second moments, namely 
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Definition 2.1 We define the Wasserstein distance between /i and u in V2 by 



W2(/i, z/) := inf 




) 
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vr : a coupling of /x and u 



(2.1) 



The function 11^2 is indeed a distance function on V2- The metric space (V2, W2) is called 
the Wasserstein space, where any two points are joined by a geodesic. Here a geodesic is 
a distance-minimizing curve and has a constant speed. It is known that any Wasserstein 
geodesic does not branch, that is, if two geodesies have a common interval, then they are 
subintervals of one geodesic. Moreover, a coupling achieving the infimum of (12. ip always 
exists, which is called an optimal coupling. For these facts, see [9[ Chapters 6,7]. We refer 
to the results of Knott-Smith ^ and Brenier [2] about optimal couplings. We denote by 
id the identity map on W^. 

Theorem 2.2 ([8, Theorems 2.12, 2.16]) Let ix^v E V2 and fi be absolutely continuous 
with respect to the Lebesgue measure. If a proper lower semi- continuous convex function 
(f) on M'^ satisfies [V0]jj/i = v, then [id xV0]j/i is a unique optimal coupling of fi and v, 
and {[(1 — t) id +tV0]jj/x}jg[o,i] is a unique Wasserstein geodesic from /i to v . 

2.2 (/^-exponential function 

We summarize definitions and properties of the 99-exponential function. For further de- 
tails, we refer to |1] and references therein. 

For an increasing, positive, continuous function ip on (0, 00), the -logarithmic function 
is defined by 



which is increasing, concave and on (0, 00). The constants l^p and L<^ are respectively 
defined as the infimum and the supremum of ln<^, that is, 

:= inf ln<^(t) = limln<^(t) G [— oo,0), := supln<^(t) = limln<^(t) G (0, +00]. 

The function ln<^ has the inverse function, which is called the ^p- exponential function and 
is defined on (Z^, L^). This inverse function can be extended to all of M as 





for r < 

In^^(T) for r G (/^,L^) , 



00 for T > L, 



which is C 



on (Z<^, L^) and its derivative is given by 



exp^(r) = (^(exp (r)). 



We mention that if p{0) = then the function exp is always 



on (-00, L^). 
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The case of ip{s) = s is the most fundamental case related to the Boltzmann-Gibbs 
statistics and the Fokker-Planck equation, where the (^-exponential function coincides 
with the usual exponential function. Another important case is that (p{s) — s^, where the 
(/^-exponential function is the power function given by 

exp^(r) := [1 + (1 - 

where we set [r]_|_ := max{r, 0} for r G M and C := oo for a < 0. This case is related to 
the Tsallis statistics and the nonlinear evolution equation of porous medium type. 
We introduce a class of increasing, positive, continuous functions on (0, oo). 

Definition 2.3 For any a G M, we define 0{a) to be the set of all increasing, positive, 
continuous functions (f on (0, oo) such that max{5^, (5'''} < a, where we set 

.1+6 

lim — -— exists /> , S'^ :— inf 5 G 



:— inf 5 G 



lim — -— = oo 



It is trivial that 0{a) G 0{b) ii a < b. We define the constant S'^, expressing the order of 
ln<^ at 0, by 



5; := inf <! (5 G 
We refer to the relation between 5^ and 5'^. 



liuit^ InJt) exists 
40 ^ 



Lemma 2.4 For any (p G 0{a) with some a G M, we have 6'^ < max{5^, 0}. 
Proof. Given any £ > 0, there exists C > such that 

„l+mEix{(5<^ ,0}+£ 

holds for any s G (0, 1). Set 6 :— max{(5^, 0} -|- 2e. Then for t G (0, 1), we have 

< -t^ InJt) = t^ f -^ds < f -^ds <C f s-^+'ds ^-(l-f)< - < +00, 
Jt Jt Jt e e 

proving 5'^ < max{5^, 0} + 2e. Since £ > is arbitrary, we have 5'^ < max{5<^, 0}. □ 



3 Condition of (p 



This section is devoted to the study of (^^-exponential distributions with mean and covari- 
ance matrix parameters. We first observe behaviors of the following function /<^. 

Lemma 3.1 For any (p G 0{1), we set the function and the constant as 



exp (A) 



(A-l%(t)r- 



t 



dt 



exp (A) 



:= max{5^, 
00 



1 if max{5^, 5^} > 0, 
otherwise. 



{X-ln^{t)y—-dt, 
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(1) For any {p, \) G (— 1,^;^) x {l^,L^), A) is well-defined. 

(2) The function A i-> f^{p,\) on {l^,L^) is C° for p G (—1,0) and for p G [0,^;^). 

(3) The function p f^{p, A) on (— l,p<^) is C°° for A G (/(^, L<^). 

(4) For any p G limA4,i^ /<p(p, A) = anc? liniAtL^ A) = oo hold. 

Proof. Note that for any p G (— l,p<^), we have l/{p + 1) > max{5^, 6'^, 6'^}. 

([1]): Fix any A G (/ip, L^) and take e G (0, A — l^). Given any p G (—1, 0), we deform 

Up^ A) = / (A - ln^(t))P--dt + / (A - ln^it)Y—-dt 

Jo fit) iexp^(A-e) V>it) 

/•exp^(A-e) ^ /.£ 

= / {X - ln^{t)f —-dt + / sPexp^(A-s)rfs 

Jo V^W Jo 

/•exp^(A-e) . 

where we use the change of variables formula for s = A — ln<^(t) in the second line, and the 
monotonicity of In^ and exp^ in the inequality. The integrability of the last line follows 
from the conditions that <f G 0(1) and p > —1. 
For any p G [0,p^), we fix 6 as 



max{5^, 5^ 5'^} <5:=l (^max{5^, 5^, 5'^} + < ^ ^ 

and deform /^(p. A) as 



/•exp^(A) ,l+<5 

f^ip, A) = {t'X - t' \n^it)Y ■ t-'^^+'^dt. 

Since 6{p + 1) < 1 holds and the function t i— )■ (t^X — t^ \n^(t)yt^^^ / (f(t) is bounded on 
(0,exp^(A)), f^{p,X) is well-defined. 

([2]), ([3]): The assertions follow from Lebesgue's dominated convergence theorem. 

(jl]): The assertion limAj,;^ fip{p, A) = also follows from Lebesgue's dominated conver- 
gence theorem. Integrating by part with the condition l/(p+ 1) > 5^ yields that 



/.(P,A) = -^(A-ln,(t))^+^ 

P -r i- 

> ^ r"p^(^) /exp^(A) - 

W(exp^(A))/ 



1 /-exp^iAj 

+ — T / (A-ln,(t))^+Mt 



P+Uo 



p+1 



_ exp^(A) / exp^(A) 
~ (p + l)(p + 2) U(exp^(A)) 

where the inequality follows from the concavity of In^, that is, 

1 m^i f t-exp (A) 
In^ t < In^ exp A + — — — = A + — -— , 

^ <^(exp^(A)) (^(exp^(A)) 
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The condition l/{p + 1) > 6"^ leads that 



exp (A) \^^^ /s^^^A^"^^ 
hm exp (A) ( — — J = hm 



AtLv \V9(exp^(A)) / sfoo I ip(s) 

proving hniAfL^ f^{p, X) = oo. □ 

Let S((i, M)+ be the set of symmetric positive definite matrices of size d, and {ei}f^i 
be the standard basis on W^. Given any V G S(c^, let V^^'^ be the symmetric positive 
definite matrix such that V^^'^ ■ V^l'^ = V and \x\y := (x, V~^x) for x E M.'^. 

Proposition 3.2 For any ip G 0{2/{d + 2)) with d > 2, there exist constants \ G {l^,L^) 
and c > such that 



exp^(A — c\x — v\y)dC = 1, (3.1) 
(x, Ci) exp^(A - c\x - v\y)dC'^ = {v, cj) , (3.2) 
r i. -v, e,) (X - „. e,) exp,(A - c|x - v\l)iC' = (e.. Ve,> . (3.3) 

Proo/. Note that the condition ip G 0{2/{d + 2)) with ci > 2 implies that 1 < d/2 < p^. 
Let r(-) be the Gamma function. Applying the change of variables formula first for 
y = c^l'^V~^l'^[x — v) then from Cartesian to polar coordinates and for t = exp^(A — |i/P), 
we have 



d-2 \ /7r\l Vdetl^ 



cJ r(c//2)' 



J ^ exp^(A - c\x - v\y)dC'^ = (—^^ ^) ' ( 

(x, Ci) exp^(A - c\x - vlDdC^ = ^) ' (^) ' ' 



{x - V, Ci) {x - V, Cj) exp^(A - c\x - vl^dC^ = (^^, A^ ■ (^^) ' ' "^^^1^^' 

which are well-defined by Lemma [3. 1[ This means that the pair (A,c) satisfying 

^((d-2)/2,A)'^ = M A-|£W2) ^ U id/2, A) 

/,(rf/2,A)i Vd^' ci/,((d-2)/2,A) ^'-^ 

is the desired one. If the continuous function F^(p, A) := f^{p — 1, \y~^^/ f^{p, Xy on 
(0,p^) X {l^,L^) satisfies 

limF^ (% X] = 0, lim (% x] = oo, (3.5) 



then the intermediate value theorem guarantees the existence of such a pair (A, c). 
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The rest is to prove (13. 5p . For any (p, A) G {—l,p^) x (l^,L^), Holder's inequality 
/•ields that 

^/,(p,A) = y^ ln(A-K(t))^-^dt 



.(P^n,A-ln,,)(^.)^(^^^^^ 

/ ^2 \ ^ 1 

providing the convexity of In in p. Then for any {p, A) G (0,p^) x (/^, L^), we have 
d d d 

— \iaF^{p, A) = (p+ 1)— ln/^(p- 1, A) +ln/^(p- 1, A) -p—\nf^{p, A) -ln/^(p, A) 

< p— In /^(p - 1, A) - p— In f^ip, A) < 
and deduce the monotonicity of in p. Since we have 

—UiO, A) = exp^(A), — /^(p, A) = P (A - ^(t))^"^^^^ = p^(p - 1, A) 

for p G de I'Hopital's rule with Lemma IXTM ]) yields that 

lim (1, A) = lim 2 exp^(A) = 0, lim F^^ (1, A) = lim 2 expj^(A) = oo, 

L p^ — \ i L p^ — \ 1 

limF,^ (p, X)p = — limF^(p - 1, A)p, lim F^ (p, A)^ = — lim F^{p - 1, A)p 

for any p G (l,p^), in particular, 

lim F^ (n, A) = lim F^ (1, A) = 0, lim F^ (n, A) = lim F^ (1, A) = cx) 

for any integer n G [l,p<p). This with the monotonicity of F^ in p finishes the proof. □ 

Remark 3.3 For the pair of constants (A, c) obtained for ^ = in (13. 4p . we find that 
exp^(A — c\x — v\y) {detV)~^^'^ becomes the probability measure with mean v and co- 
variance matrix V by the change of variables formula for y = V^^l'^ix — v). 



4 Convexity in Wasserstein geometry 

In this section, we investigate the two spaces of (^-exponential distributions with mean 
and covariance matrix parameters. We mention that a^pi^s) = s'' G (9(2/((i + 2)) with 
some a > if and only if g G := (0, {d + 4)/(ci + 2)). 



7 



Proposition 13.21 and Remark 13.31 yield that for any ip G (9(2/(c/ + 2)) with d>2, there 
exist continuous functions and on S('i, such that each of 

n^{v, V){x) := exp^ (A^(V) - c^{V)\x - v\y) , 

g^{v, V){x) := exp^(A^(/rf) - c^{Id)\x - v\y){detV)-^ 

is the probabihty density on {W^, C^) with mean v and covariance matrix V. We define 
the two spaces of probabihty measures by 

Af^ := {N^iv, V) := n^{v, V)^ \ (v, V) eR'' x SK 

:= {G^{v, V) := g^iv, V)C'' \ (v, V) eR'' x S(rf, M)+}, 

which generally differ from each other. 

Theorem 4.1 For any <f,ipE 0{2/{d + 2)) with d>2, G^{0,a^ld) = N^{0,a^ld) holds 
for any a > if and only if aip{s) = f3ip{s) = s'^ with some a, /3 > and q E Qd- 

Proof. We first prove the "if part. For any Lp E 0{2/{d + 2)) with d>2 and a > 0, the 
following relations are verified; 

ln„^(t) = ln<^(t), exp^^(r) = exp^(ar), A) = a"^^+^V</p(p, "A). 

These relations with (13. 4p implies that, for any V E S('i, we have 

W-f^^ = F^, {^A.,{V)^ = F, {^,aK,{V) 

fa^{d/2,K^{V)) f^{d/2,aX^^{V)) 

Coup \' j — , — OL 



df^^ {{d - 2)/2, Xa^V)) df^ {{d - 2)/2, aX^^iV)) ' 

which shows that aXaipiV) = X^p{V) and aCa^iV) = c^p{V), in turn Ga,fi{v, V) = G^{v, V) 
and Naip{v,V) = N^{v,V). Therefore it is enough to prove the case (p{s) = ipi^s) = s'^ 
with q E Qd, which has been already demonstrated in [71 Section 4]. 

Let us prove the "only if part. The condition G^(0,a^/rf) = N^{0,a'^ld) leads that 

a~'^exp^(A — ca~^|x|^) = exp^(A(a) — c(a)a^^|x|^) 

for any x E M.'^, where we abbreviate c := c^{Id), X := X^{Id) and c(a) := c^{a^Id), A(a) := 
X^[a^Id)- Evaluating the above equation at a; G M'' with = ar, we deduce that 



a-'^exp,(A - cr^) = exp^(A(a) - c(a)r2), h := = ^ 

Differentiating this equation at r G / := (0, h) yields that 

ca~'^V9(exp^(A — cr^)) = c(a)?/^(exp^(A(a) — c(a)r^)) = c(a)'?/'(a~'^exp^(A — cr^)), 
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in particular, ap{exp^{\ — cr^)) = c(l)'?/'(exp^(A — cr^)) for the case a = 1. This provides 

V>(a-'^exp^(A-cr^)) _ c(l)^(a-^ exp^(A - cr^)) _ c(l) ^ 
'?/'(exp^(A — cr^)) cv9(exp;^(A — cr^)) c{a)a'^ 

for any r E I. For any ^,r7 > satisfying ^^,^?7 < expj^(A), since it is possible to choose 
s,r E I such that = exp^(A — cr^) and = exp^(A — cs^), the choice = v/^ leads 

V^(e^) ^ V^(a"^exp^(A-cr=^)) ^ c(l) ^ V'lr^^) 
^(^2) ^(exp^(A - cr2)) c(a)a°' V'(^?7)' 

Then for any a E (0, min{expj^(A)"'^/2, expj^(A)}), the function '■= 
satisfies that ^{^r]) = \E'(^)\I'(?7) for any ^,r] E {0, a"^^"^ exp^^Xy^"^) . By induction with 
the condition a < exp^(A), we find that ^^(^r]) = \E'(^)\E'(?7) for any ^,77 > 0. Hence the 
continuous function L{s) := ln\E'(e'*) on M satisfies Cauchy's functional equation, namely 
L{s + 1) = L{s) + L(t). This means that L{s) = L{l)s, in turn aip{s) = f3ip{s) = s'' for 
some a, (3 > and q E Qd- n 

Remark 4.2 Theorem 14.11 yields that for ip^ip E 0{2/{d + 2)) with d > 2, we have 
^^(0, a^Id) = N^{0, a? Id) for any a > if and only if and are either the space of 
Gaussian measures, or the space of g-Gaussian measures. 

Both the spaces Q^p and A/"^ are generalizations of the space of Gaussian measures, however 
they differ from each other in general. On one hand, the space has a dually flat 
structure as seen in [4] (we do not refer to a dually fiat structure, see [Tj). On the other 
hand, the space is convex in Wasserstein geometry. 

Proposition 4.3 For any ip E 0{2/{d + 2)) with d > 2, the space is convex and 
isometric to the space Q of Gaussian measures in the sense of Wasserstein geometry. 

Proof. For (f , V), {u, U) E M'^x S((i, + , we define the symmetric positive definite matrix 
W and the associated convex function cpw by 

W := U^iuWu^^Y^^U^, (t)w{.x) := ^{x - v,W{x - v)) + {x,u). 

Then (p^y is convex and satisfies [V4>w]iGy,{v,V) = G^{u,U) since for G^p{v,V)-a\most 
everywhere we have g^p{v,V) = g^{u, U) (V (pw) {det B.ess (pw)- Theorem 12.21 vields that 
[id xV(l)w]iG^{v, V) is an optimal coupling of Gif,{v, V) and G^{u, U), consequently 

W2{G^{v,V),G^{u,U)y = |i;-M|2 + trl^ + trf/-2tr (u^VU^Y 

= W2{G{v,V),G{u,U)y, 

where G{v, V) stands for the Gaussian measure with mean v and covariance matrix V. 
(Note that the optimality of the function (pw in the case of the Gaussian measures has been 
known, for instance, see [6] and references therein.) Hence the map G{v, V) 1— )■ G^{v, V) 
is an isometry from Q to Q^p in Wasserstein geometry. We also deduce from Theorem 12.21 
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that, for the time-dependent vector {wt}te[Q,i] and the time-dependent matrix {iyj}jg[o,i] 
defined by 

wt := (1 - t)v + tu, Wt := [(1 - t)/^ + tW]V[{l - t)h + tW], 

{G^p{wt^Wt)}t(^[Q^i] is a unique Wasserstein geodesic from G^iy^V) to G^{u,U). Hence 
the space is convex in Wasserstein geometry. □ 

Let us now consider the completions of as a metric space with respect to W2, 
denoted by 7\7^, and prove that J\f^ is generally not convex in Wasserstein geometry. By 
Proposition 14. 3[ the completion of Q^, denoted by Q^, is isometric to the completion of 
Q, which is homeomorphic to M.'^ x S('i, K)>0) where S('^, K)>o is the set of all symmetric 
no n- negative definite matrices of size d (see ^ Section 4]). We similarly show that A/^ is 
homeomorphic to M'^xS((i, K)>0; whose proof is simple but tedious and we omit it. The key 
of the proof is the fact that the weak convergence, which is equivalent to the pointwise 
convergence of characteristic functions, is weaker than the convergence in Wasserstein 
geometry (see [8l Theorem 7.12]). Hence elements of and A/^ are naturally denoted 
by G^{v, V) and N^{v, V) with some {v,V) eW^ x S{d, M)>o, respectively. Note that any 
elements in \ and \ A/"^ are singular with respect to the Lebesgue measure and 
G^{v, 0) = N^{v, 0) = 6v, which is the Dirac measure centered at v. 

Proposition 4.4 The space is convex in Wasserstein geometry if and only ifaip{s) = 
s'^ with some a > and q ^ Qd- 

Proof. Since the "if part is trivially true, we assume the convexity of A/^. Then a unique 
geodesic a from A^^(0, 0) = G^(0, 0) to iV^(0, Id) = G^{0, Id) is contained in A/^ and Q^, 
which is written as a{t) = G^{0,t'^ld) for t e [0, 1] and can be extended to [0, 00). The 
non-branching property in Wasserstein geometry ensures that G^(0,t^/d) G fl A/"^ for 
any t > and by Theorem 14.11 aip{s) = s"^ holds for some a > and q G Qd- ^ 



5 Stability under an evolution equation 

We finally discuss when the space of y^-exponential distributions is stable under the evo- 
lution equation of the form 

= div (pV(ln^(p) + ^^)) = div + pV^^^ , ^^{x) = c^{Id)\x\^. (5.1) 

We mention that this evolution equation recovers the Fokker-Planck equation (resp. the 
nonlinear evolution equation of porous medium type) when (f{s) = s (resp. if{s) = s'^ 
with some q & Qd except for g = 1). 

A set of functions is said to be stable under an evolution equation if any solution to 
the evolution equation with initial data being an element of the set stays inside the set. 
A set of absolutely continuous measures with respect to the Lebesgue measure is said to 
be stable under an evolution equation if the set generated by densities of elements in the 
given set is stable. 
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Theorem 5.1 Assume y:>{0) = 0. Then the space M^p is stable under (15. ip if and only if 
aip{s) = s'' with some a > and q G Qd- 

Proof. Since the "if part has been aheady shown in [21 Proposition 5], it is enough to 
prove the "only if part, which is done in a similar way to the proof of [S]. 

Assume the stability of Af^. Then the solution to (15.11) with initial data n<^(0, V) G J\f^ 
is written as n^{vt, Vt) for some time- dependent vector vt G M'^ and some time-dependent 
matrix Vt G S(c?, The following calculation can be justified; 



d f 
^T^n^i^t, VtjdC^ = - n^ivt, Vt)V (ln^(n^(ft, Vt)) + dC^ = -2cp{Id)vt, 



d 

dt""^- 

which implies that Vt is identically the zero vector. We similarly compute 

jVt = AAtVt, jVt-' = -Wt-'At, At := c{Vt)Vt-' - Cp{Id)U (5.2) 

In what follows, we abbreviate n^{0,X) as n^p{X). Let dij be the tangent vector on 
J\fip{0) := {n^{X) I X G S(rf, IR) + } associated to a global chart (S((i, M)+,S = {'^ij)i<i<j<d) 
given by S(X) := n^(X~^). We define the map a on S('i, K)+ by cr(X) = X~^ and the 
symmetric matrices A, C of size d by 

{<9ij(A<^ocr)|^_i iii=j, ( dij{cp o a)\^-i if ^ = J, 

1 ! C := \i ' 

Given any x G M*^, we set the function Cx on A/'^(0) as ex{ny,{X)) = np{X){x), which is 
owing to f{0) = 0. Then for the solution np{Vt) to (15.11) . we have 

^e^MV)) = ^npiVt)ix) = div (n^(\4)V(ln^(n^(V;)) + ^^)) (x) 

= Ay^{np{Vt){x))c{Vt) {x,Vf'Atx) - 2n^(\4)(x) tr(A). 
On the other hand, by using (15. 2p . a direct computation provides 
d , ^ ^ I - \\ d{ei,Vt-^ej) 

= -Aip{np{yt){x)) [tiiYt-^Atk) - ti{yt-^AtC)\x\\ - Cpiyt) (x, Vt-'Atx)] . 

We therefore obtain 

np{Vt){x)tT{At) =2ip{np{Vt){x)) [tr(Vr'AA) - tr(Vr'AC)|x|2.)] (5.3) 

for any a; G M'^. Set 

_ 2tTiVt-'AtA) . 2tT{Vt-'AtC) , . X nn l \ 

^ ■= triAt) ' c,(T4)tr(A)' Jv.-=(^^KiV.)-l,). 
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At X G M'^ with C(^(Vf)|x|y^ = r G Jvj, the equation (15. 3p is deformed as 
exp^(A^(Vi) - r) = y.(exp^(A^(ri) - r)) (/3 + 7r) = - f exp^(A^(Vi) - r) ) (/3 + ^r) 



and integrating it with respect to r G Jvt yields 



exp^(A^(Vi) - r) = exp^(A^(V-0) (l + 



Taking In^ and differentiating it in r, we have 

/3exp^(A^(rt))Xs) = si+^ 

where we set s = s{r) := exp^(A^(Vt) — r). This imphes that aip{s) = s'^ holds on 
(0,exp^(A(\/i))) with a = a{Vt) := (3 exp^{X^{Vt))"' and q = q{Vt) := 1 + 7- We mention 
that a and q do not depend on Vt since we have (p{s) = a{Vt)~^ s'^^^''^ = a{Ut)~'^ s'^^^*^ 
on (0,min{exp^(A(\4)),exp^(A(f/t))}) ^ for any V,U e S(rf,M)+. Letting V ^ and 
t — > 0, we have A^(Vt) — )■ then 0^9(5) = s'' for any s > 0. □ 

The same result holds for Q^,. The proof is exactly similar but tedious and we omit it. 



Corollary 5.2 Let ip be and (p{Q) =0. Then the space Q^p is stable under (15.11) if and 
only if aip{s) = with some a > and q G Qd- 
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